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We study the continuum field theory for an ensemble of directed polymers 
Ti(t) in 1 + d' dimensions that live in a medium with quenched point disorder 



?-j | and interact via short-ranged pair forces ^(rj — r^). In the strong-disorder 

(or low-temperature) regime, such forces are found to be relevant in any 
dimension d! below the upper critical dimension for a single line. Attractive 
forces generate a bound state with localization length £j_ ~ Isl - "- 1 ; repulsive 
forces lead to mutual avoidance with a pair distribution function V{ti —fj) ~ 
|rj — rj\ s reminiscent of interacting fermions. In the experimentally important 
dimension d! = 2, we obtain z^j_ ~ 0.8 and 8 ~ 2.4. 
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Flux lines in dirty type-II superconductors |1[] are a well-known example of low- 
dimensional manifolds embedded in media with quenched disorder (reviewed e.g. 
in p|). Such systems are of interest also because of their links to more complicated 
random systems such as spin glasses ||, to surface growth |4], [J, and to randomly 
driven hydrodynamics ||. A well-studied case is the low-density limit of a single 
line in a medium with quenched point impurities. The statistical properties of such 
a line differ from those of a free, thermally fluctuating line. Moreover, the disorder 
modifies the interactions of the line with other objects. For example, the effect of a 
rigid line defect ("columnar" defect) parallel to the preferred axis of the lines turns 
out to be weaker than in a pure system. A weakly attractive defect localizes the line 
only up to the borderline dimension d'^ — 1; in higher dimensions, the transition to 
a localized state takes place at finite coupling strength |7|, || . In a pure system, the 
borderline dimension is d^ = 2 ||. 

Mutual interactions between several identical lines in a disordered medium are 
the subject of this letter. We study their effects by the methods of continuum field 
theory. Related aspects of this system have been treated by Bethe ansatz methods 
in d' — 1 [|ITJ, on a lattice [O, O, and in a Wilson renormalization group [|TJ, KM , 
see the discussion below. 

In a pure system, weakly attractive short-ranged forces generate a bound state 
up to the same dimension d'+ = 2 as do columnar defects. Unlike in the case of 
columnar defects, however, the effect of pair interactions between identical lines 
is enhanced by the disorder, which (at sufficiently low temperature) leads to a 
bound state in any dimension where the low-temperature behavior of a single line 
is governed by non-thermal scaling exponents. Of equal importance are repulsive 
forces; for example, the magnetic interaction between flux lines can be treated as 
a contact interaction at sufficiently low densities |IJ, [Tj|. In a pure system, such 
forces are important only in d! = 1, where they act as an effective constraint on 
the fluctuations that is equivalent to the Pauli principle: the lines behave like the 
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worldlines of free fermions characterized by a pair distribution function P(rj — Tj) ~ 
\t{ — Tj\ e with 6 = 2 |y| 16| . In higher dimensions, there is no long-ranged effect on 



the pair-distribution function. In the presence of disorder, on the other hand, we 
find that Viji — Tj) obeys a power law also for d' > 1, with a new exponent 6 that 
depends on the dimension. This behavior is reminiscent of interacting fermions, 
which is why we call these mutually avoiding lines fermions in a random medium. 
This strong effect of repulsive interactions can be understood qualitatively as 
resulting from the energetic competition with an effective line-line attraction due to 
the impurities. In almost all realizations of the disorder, there is a unique ground 
state, i.e. a path of minimal energy ||17|| . At low temperatures and without direct 
forces, two lines will share this path with finite probability even in a system of 



infinite size pA\ , [19], P0|| . Repulsive interactions, however, force one of the lines into 
a (distant) excited state. 

Our calculations in the sequel will be restricted to the case of two interacting 
lines in a random medium; the generalization to an arbitrary number of lines is 
straightforward. On a mesoscopic scale, the system is described by an effective 
continuum Hamiltonian 
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Here ri(£) and ^(i) denote the displacement vectors of the two lines (also called 
directed polymers) in d! transversal dimensions, as a function of the longitudinal 
"timelike" coordinate t. Both lines are subject to the same random potential 7/(r, t) 
which models the quenched point defects. It is assumed to be Gauss-distributed 



with r](r, t) = and rj(r, t)r](r', t') = 2<r 2 S d (r — r')5(t — t'). Averages over disorder 
are denoted by an overbar, thermal averages by brackets (•••). The direct inter- 
action between the lines, which is assumed to be short-ranged, is described by the 
continuum field ^(t) = 5 d ' (r x (t) - r 2 (t)). 

For g = 0, the two lines feel only the random potential, and the free energy 



of the system is just twice the free energy of a single line. The continuum field 
theory for this system serves as the starting point for a systematic renormalized 
perturbation theory in the line-line coupling g. The analysis is complicated by the 
fact that even the "free" theory (g = 0) has non-Gaussian multipoint correlations 
due to the quenched averaging. Nevertheless, we are able to obtain solely in terms 
of single-line properties the scaling dimension and the form of the operator product 
expansion of the pair interaction field *&(t), see (J3|) and (U) below. These determine 
the renormalization group equations for the interaction strength to leading order, 



and hence the phase diagram of the system [21]. 

The large-scale behavior of a single line is generated by the sample-to-sample 
fluctuations of the ground state paths and defines two important exponents. The 
roughness exponent ( characterizes the mean transversal excursions of the line, given 



e.g. by the two-point function ((r(t) — r(£')) 2 ) ~ \t — t'\ 2 ^. The exponent — u is the 



anomalous dimension of the disorder-averaged free energy F = /3 _1 logTrexp(— /3TC), 
whose universal part has the scaling form F(T, R) ~ T w jF(i?/T < ') in a finite system 
of transversal size R and longitudinal size T |[22|| . The two exponents are related by 
a "tilt" symmetry of the system, u = 2( — 1 (see e.g. [IS, Kfl). In low dimensions, 
( is always larger than in the case of thermal fluctuations, namely ( = 2/3 for 
d' = 1 and £ « 0.62 for d! = 2. For d' > 2, a phase transition appears at a 
finite temperature; in the high-temperature phase, the system is asymptotically 
thermal, i.e. ( = 1/2 and uj = (see e.g. || and references therein). Whether a 
finite upper critical dimension d> exists such that for d' > d^ the thermal exponents 
govern also the low-temperature phase is controversial; some workers believe d> ~ 4 
2l|. In the continuum theory (JJ) with g = 0, the large-scale regime is reached 



in a crossover from the Gaussian theory with characteristic longitudinal length 
£u = P(a 2 p 3 )^ 2 ^ 2 ^ d \ We have discussed the renormalized continuum field theory 
for this regime in ref. ||. Its construction involves a reparametrization of the 
transversal displacement and of the free energy, r — > P 1 ^ 2 ^ r and F — > /3&F, 



such that the renormalized variables remain finite in the limit £u — » (i.e. (3^ 1 — > 
or cr 2 — > oo). The additional reparametrization \1/ — > P~ d '^ 2 t\\ ^ an d g — * 
^l+d /2£w ^ keeps all correlation functions of the pair field free of singularities 
in this limit, as we will see explicitly below. This is accompanied by a temperature 
reparametrization; the renormalized temperature f3~ l ~ & is an irrelevant scaling 
variable of dimension —u. 

With this renormalized theory for g = at hand, the effect of pair interactions 
may now be expanded in powers of the coupling constant g. For example, the 
perturbation series for the free energy density / = limr^oc, F/T in a system of 
transversal size R = Z>, 

:-Pg) N 

JV=1 

involves integrals over connected correlation functions of the pair field \l/ evaluated 
at g = 0. We are hence led to study these correlation functions. While they are 
finite by construction, the "time" integrations generate new singularities that have 
to be handled by an additional renormalization. (The longitudinal scale L plays 
the role of an infrared cutoff and will be used to parametrize the renormalization 
group flow.) 
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The one-point function (&(£)) gives the probability density that the two lines 
intersect at time t, averaged over thermal and disorder fluctuations. We will show 
below that for R 1 ^ > £,, , 



mt))=bR°, (3) 

where b is a constant independent of £n. Hence \l/(t) can be regarded as a scaling field 
of dimension zero. Eq. (^) expresses the fact that for g = 0, the two lines follow the 
same path (of minimal energy) with finite probability even in the thermodynamic 
limit R — >• oo. 

The higher connected correlation functions of the local pair field \1/ can be shown 



to obey the operator algebra 

V(t)*(f) = cfi- l \t - t'^ V(tf) + ... (4) 

with a coefficient c > 0, which is valid as an asymptotic identity inserted into 



any such correlation function (. . . \l/(t)\l/(t') . . .) c for t — > if. This type of operator 
algebra is familiar from ref. flSJ (where the reader is referred to for a more detailed 
discussion): the field \1/ couples to itself, but the leading singularity in (f|) involves 
a correction-to-scaling exponent related to the irrelevant coupling constant (3~ x . 
To establish eq. (|3]), it is useful to study the pair distribution function V(d) = 



(S d '(ri(t) — r 2 (i) — d)), i.e. the probability density that the two lines have the 
relative displacement d. We find that for g = and in the thermodynamic limit 
R — > oo, V has a scaling form that depends only on |d| and £y with the asymptotic 

behavior 

^| d |-KC+-)/C for |d|»| l| C , 

Z\r d ' ( for |d|<^„ , 

see also the extensive discussion in ||20|| . Recall that ^ is the renormalized tem- 
perature /9 _1 . Hence with finite probability, the two lines share a common "tube" 
of width £[[ (i.e. |d| < Q), but they do make large excursions whose probability 
decays with a power of |d|. To show this, we first rewrite V(d) for g = and d 7^ 
in terms of the properties of a single line, 

V(d) = fd d 's (5d'( ri (t)-s)5 d '(r 2 (t)-s-d)) (6) 

= Jd d 's(<5>(t,s)) l ($(t,s + d)>! (7) 

= -|d d 's($(t,s)$(t,s + d))f . (8) 

In the single-line system, $(£, s) = 5 d (r(t) — s), and averages are marked by the 
subscript 1. Here we have exploited the fact that without pair forces, the two lines 
are independent in any realization of the disorder, and that the full correlation func- 
tion ($(£, s)$(£, s + d))i vanishes since a directed line has a single- valued path r(t). 



Next we use the operator algebra $(£, s)<3>(£, s + d) ~ -ci/? _1 |d| ( rf/ M-«0/C$(£ ; s ) + . . . 
with c\ > (that follows directly from the results of |§ or [^Oj) and the normal- 
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ization / d d s (<&(£, s))i = 1 to obtain ([5]) in the regime |d| ^> £y ; the oppo- 
site limit then follows from the normalization / d d d V(d) = 1 and the relation 
Jd d d V(d) d 2 = 2f3~ 1 R 1 ^, a consequence of the tilt symmetry |24j| . For R — > oo, 



notice that P(d) as given by (|5]) is only normalizable for uj > 0. The expectation 
value of the pair potential \I/ at g = involves an integral over the potential func- 
tion times V(d). For a potential of microscopic range Ao which is well-behaved at 
d = 0, this integral is independent of R, which is just eq. (|3|). For A <C Q, the 
integral develops a singularity £y ^ that we have absorbed into the definition of the 
renormalized field \I/. 

The form of the operator algebra (§) could be derived by arguments very similar 
to the ones we have used in ref. ||. Here we take a different avenue: we exploit 
the exact mapping of the quenched polymer system defined by ([!]) to a nonlinear 
stochastic evolution equation, a generalization of the Kardar-Parisi- Zhang equa- 
tion |jj. If Z(ri,T2,t) denotes the restricted partition sum over all two-line config- 
urations with fixed end points (ri,r 2 ,£), the field h(r 1 ,r 2 ,t) = j3~ 1 \ogZ(ri,r 2 ,t) 
obeys the equation 

dth = u £V 2 ^ + ~ E(V r ^) 2 +J2r,(r i ,t)-g8 d '(r 1 -r 2 ) . (9) 



Mukherji [|l l| uses this equation to treat the problem of two interacting lines by 
a standard dynamic renormalization group approach, expanding simultaneously in 
g and A (see the further remarks below). From (|9|), one constructs in a stan- 



dard way the generating functional of the dynamic correlation functions |25|| , which 
are denoted by ((•••)). This functional integral involves the field h and an aux- 
iliary field h. Insertions of the auxiliary field generate response functions, e.g. 
((/i(ri,r 2 ,t)n™ = ifc(rifc,r 2fc ,t fc )>) = S m ((h(r u r 2 ,t)))/UT=i^J^ik,r 2k ,t k ), where 
J(ri,r 2 ,t) is a source term added to the r.h.s. of (§). By functional differenti- 



ation with respect to J, one easily verifies the following relation between connected 
averages in the polymer system and dynamic response functions, 

/T" 1 WWh) - si)^'(r 2 (ii) - si) . . . <K(n(£ m ) - s m )<K(r 2 (£ m ) - OF = 

lim {(h(r 1 ,r 2 ,T)h(s 1 ,s' 1 ,t 1 )...h(s m ,s' m ,t m ))) . (10) 

T — >oo 

The dynamic average ((h(r 1 , r 2 , t))) equals (minus) the disorder-averaged free energy 
of the polymer system. For g — 0, one has ((/i(ri, r 2 , t))) = Ya=i ((M r J'*)))ii hence 
/i(r l5 r2,t) is a scaling field of dimension — uj (with time as the basic scale), as in 



the usual Kardar-Parisi-Zhang equation. From (|10D, we have in particular (^(t)} = 
HmT^oo((h(ri,r2,T) ip{t))) with ip(t) = Jd d r h(r,r,t). Comparing dimensions on 
both sides, we conclude that ip(t) is a dynamic field of dimension uj. Hence its 
self-coupling has to be of the form ij;(t)ip(t') ~ c'\t — t'\~ u ip(t) + . . .. Transforming 
back to the polymer system then gives @. 

The renormalization of the perturbation series (0) to one-loop order is based 
on eqns. (^) and (|j); it can be borrowed entirely from ref. M. It yields the flow 
equation LBlUr = eur — cu 2 R + 0(u R ) for the renormalized pair coupling ur (that 
is proportional to g), with e = 1 — oj(d'). To the given order, this equation has the 
two fixed points Ur = and u* R = (1 — u)/c. In any physical dimension d', one 
has e > 0. Hence the trivial fixed point ur = is unstable, and the interaction is a 



relevant perturbation. This result is at variance |g6| with ref. [13 . 

Any attractive line-line force (g < 0) grows indefinitely in magnitude under the 
renormalization. The system develops a bound state with transversal localization 
length £j_. As g approaches zero from below, the lines unbind continuously, i.e. the 
localization length diverges as 

t±~\g\~ u± with i/ ± = C/e=(l + w)/2(l-a;) (11) 

when the interaction strength approaches zero. 

With a repulsive interaction (g > 0), however, the large-scale behavior of the 
two-line system is determined by the nontrivial fixed point u* R . At this fixed point, 



the pair field \1/ acquires the new dimension x* = 2(1 — uo) + 0((1 — uo) 2 ) that de- 
termines the large-scale asymptotics of the pair distribution function. For example, 
the probability of intersection in a system of transversal size R scales as 



(#(£)) ~ R~° with 6 = x*/( , (12) 

where 9 « 2 for d' = 1 and 6 1 « 2.4 for d' = 2. Unlike in the case without pair 
interactions, the probability of intersection now approaches zero in the thermody- 
namic limit R — > oo: the lines avoid each other completely. Moreover, this implies 
by standard scaling arguments that repulsive pair forces have a long-ranged effect 
on the pair distribution function, 

V{v 1 -v 2 )^\v 1 -v 2 \ e (13) 

for i?> |n -r 2 | »£l(sO- 

We may use these results to obtain the finite-size scaling of the "overlap" 
q(T,g) = T _1 / dt' (ty(t'))(T, g) in a system of longitudinal length T and with 
R — > oo, a quantity that is discussed in the literature fill O, [14|. We have 



q(T,g) = q(T,0)Q(T/Q = Q{g T 1 -), using g(T, 0) ~ T° by © and the seal- 
ing of the logitudinal correlation length £u ~ £/ given by (|11~D. For d' = 1, where 



1 — uo = 2/3, this scaling form agrees with Mezard's conjecture |TT| (which is based 
on numerical simulations) and with Mukherji's renormalization group results [[14 



The latter approach, however, gives no information on the strong coupling phase in 
any dimension d' > 1. 

In summary, we have shown that directed lines in a highly disordered medium 
respond to pair forces in a stronger way than lines with purely thermal fluctuations. 
In any dimension d' < d' y , they form a bound state with attractive forces; with 
repulsive forces they avoid each other, as described by the pair distribution function 
flTB"!). The reason for this strong effect of short-ranged interactions is that when they 
are absent, the disorder forces the lines to cluster in the vicinity of a unique path 
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of minimal energy. Our results hence provide an experimentally and numerically 
accessible consequence of the clustering. However, the pair distribution function 
(H) shows a singular broadening for fixed £y and R — > oo as d! approaches the upper 
critical dimension d' y (i.e. ui \ 0): the width A(rj) defined by ij,j|<A d d d ^(d) = 77 



diverges as A (77) ~ l/u> for any < 77 < 1, and accordingly (\I/) ~ to. If d! > is finite, 
we hence expect that for d! > <i>, the lines no longer cluster, but exploit multiple 
near-minimal paths even at low temperature, as in a glassy state. This leads to a 



modification of eq. (|3|), (\I / ) ~ R~ d ; weakly attractive pair forces should then no 
longer generate a bound state. This could be useful to determine d! > numerically 
in a way that is less hampered by finite-size effects than the existing simulations of 
the KPZ equation 



References 

[1] M.V. Feigel'man, V.B. Geshkenbein, A.I. Larkin, and V.M. Vinokur, Phys. Rev. Lett. 63 
(1989), 2303; T. Hwa, Phys. Rev. Lett. 69 (1992), 1552; D. Nelson, H. S. Seung, Phys. Rev. 
B 39 (1989), 9153; D. Nelson, P. Le Doussal, Phys. Rev. B 42 (1990), 10113. 

[2] G. Forgacs, R. Lipowsky, and T.M. Nieuwenhuizen, in Phase transitions and Critical Phe- 
nomena, Vol. 14, ed. C. Domb and J. Lebowitz (Academic Press, London, 1991); M. Kardar, 
Lectures on Directed Paths in Random Media, in Fluctuating Geometries in Statistical Me- 



chanics and Field Theory, Les Houchcs proceedings, preprint cond-mat/9411022 



[3] B. Derrida and H. Spohn, J. Stat. Phys. 51 (1988), 817; M. Mezard and G. Parisi, J. Phys. 
I (France) 1 (1991), 809; H. Kinzelbach and H. Horner, J. Phys. I (France) 3 (1993), 1329 
and 1901. 

[4] M. Kardar, G. Parisi, and Y.-C. Zhang, Phys. Rev. Lett. 56 (1986), 889. 

[5] J. Krug and H. Spohn, in Solids Far From Equilibrium: Growth, morphology and Defects, 
ed. C. Godreche (Cambridge, 1990). 

[6] D.A. Huse, C.L. Henley, and D.S. Fisher, Phys. Rev. Lett. 55 (1985), 2924; D. Forster, D.R. 
Nelson, M. Stephen, Phys. Rev. A 16 (1977), 732. 

[7] M. Kardar, Nucl. Phys. B290 (1987), 582; L. Balents and M. Kardar, Phys. Rev. B49 (1994) 
13030; L.-H. Tang and I.F. Lyuksyutov, Phys. Rev. Lett. 71 (1993), 2745; T. Hwa and T. 
Nattermann, Princeton report, |cond-mat 9404031 (1994). 



[8] H. Kinzelbach, M. Lassig, Depinning in a Random Medium, preprint cond-mat/9405088 

[9] R. Lipowsky, Europhys. Lett. 15 (1991), 703. 
[10] M. Kardar and D.R. Nelson, Phys. Rev. Lett. 55 (1985), 1157. 
[11] M. Mezard, J. Phys. (France) 51 (1990), 1831. 
[12] L.-H. Tang, J. Stat. Phys. 77 (1994), 581. 

[13] T. Nattermann, M. Feigelman, I. Lyuksyutov, Z. Phys. B84 (1991), 353. 
[14] S. Mukherji, Phys. Rev. E 50 (1994), R2407. 

10 



[15] H.B. Thacker, Rev. Mod. Phys. 53 (1981), 253; M.E. Fisher, J. Stat. Phys. 34 (1984), 667. 

[16] S. Mukherji and S.M. Bhattacharjee, J. Phys. A 26 (1993), L1139; M. Lassig, Phys. Rev. 
Lett. 73 (1994), 561. 



[17] Due to the "tilt symmetry", one has ((r.j(£) — r.;(£')) 2 ) c = /3 _1 \t — t'\ exactly. A degeneracy of 
ground states in an extensive number of realizations would lead to temperature-independent 
corrections to this relation. See the discussion in [Q g(J and the results of |L2] , 

[18] D.S. Fisher and D.A. Huse, Phys. Rev. B 43 (1991), 10728. 

[19] G. Parisi, J. Phys. (France) 51 (1990), 1595. 

[20] T. Hwa, D. Fisher, Phys. Rev. B 49 (1994), 3136. 

[21] M. Lassig and R. Lipowsky, Phys. Rev. Lett. 70 (1993), 1131, discuss a similar renormaliza- 
tion group for temperature-driven transitions of directed manifolds. 

[22] D.A. Huse and C.L. Henley, Phys. Rev. Lett. 54 (1985), 2708. 



[23] See e.g. M.A. Moore et al., preprint :ond- mat/9407076, and references therein 



[24] /d d d V(d) d 2 = ((n(t)-r 2 (t)n = 2((r(t))' )f = hm,^,^ <(r(t) - r(f))' )?. For 



\t — t'\ <C R 1 ^ , this equals (3 x \t — t'\ due to the tilt invariance ITq, 20|; in the opposite 



limit, one then obtains (3 l R 1 '^ by finite-size scaling arguments; see also pOfl . 

[25] For the general formalism, see e.g. P.C. Martin, E.D. Siggia, H.A. Rose, Phys. Rev. A8 
(1973) 423; R. Bausch, H.K. Janssen, H. Wagner, Z. Phys. B24 (1976), 113. Applied to the 
KPZ-cquation, sec T. Sun, M. Plischkc, Phys. Rev. E49 (1994), 5046; E. Frey, U.W. Timber, 
Phys. Rev. E50 (1994), 1024. 

[26] We attribute this discrepancy to the fact that the disorder-induced attraction between the 
lines is neglected in the treatment of ]13[ ]. 

[27] T. Ala-Nissila et al., J. Stat. Phys 72 (1993), 207, and references therein. 



11 



